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After reviewing QCD definitions of the chiral-odd spin-dependent parton distribu- 
tions hi(x, Q 2 ) and h^(x, Q 2 ), I will summarize the main feature of the recent two 
results in perturbative QCD: (i) Next-to-leading order Q 2 evolution of hi(x, Q 2 ). 
(ii) Leading order Q 2 evolution of the twist-3 distribution h^{x,Q 2 ) and the uni- 
versal simplification of the Q 2 evolution of all the twist-3 distributions in the large 
N c limit. 



1 Introduction 

Recent developments in collider experiments have been providing us with rich 
information on the quark-gluon distributions in the nucleon. Particularly inter- 
esting are the spin structure functions which reveal "spin distributions" carried 
by quarks and gluons inside the nucleon. Besides the importance in their own 
right, they play an indispensable role to test the spin-dependent part of QCD - 
fundamental theory of the strong interaction: Deeper test of QCD in the spin 
dependent level may eventually lead us to search for a new physics beyond the 
standard model. 

The nucleon's structure functions measured in hard processes can be writ- 
ten as a sum of parton distributions for each quark or anti-quark flavor and 
a gluon. They are functions of Bjorken's x which represent parton's momen- 
tum fraction in the nucleon and a scale Q 2 at which they are measured. Un- 
till now, most data on the nucleon's structure functions have been obtained 
through the lepton-nucleon deep inelastic scattering (DIS). The chiral-odd dis- 
tributions, hi l(x, Q 2 ), however, can not be measured by the inclusive DIS, 
and hence there has been no data up to now. They can be measured by the 
nucleon-nucleon polarized Drell-Yan process and semi-inclusive DIS which de- 
tect particular hadrons in the final state. They will hopefully be measured by 
planned experiments using polarized accelerators at BNL, DESY, CERN and 
SLAC etdil. In particular, RHIC at BNL is expected to provide first data on 
these distributions. 



"Invited talk presented at "QCD Corrections and the New Physics", October 27-29, 1997, 
Hiroshima, Japan. To be published in the proceedings, (ed. by J. Kodaira et al.) 



1 



In the study of these structure functions, perturbative QCD plays an in- 
dispensable role in predicting their Q 2 -dependence: Given a structure func- 
tion, say hi(x,Qo), at one scale Qq, perturbative QCD predicts the shape of 
hi(x, Q 2 ) at an arbitrary scale Q 2 . This Q 2 evolution is necessary not only in 
extracting low energy hadron properties from high energy experimental data 
but also in testing the x-dependence predicted by a non-perturbative QCD 
technique or a model with the high energy data. In this talk, I will first give 
the QCD definition of hiL(x,Q 2 ) and discuss their Q 2 -dependence studied in 
the recent literature. 




Figure 1: (a) Quark distribution function, (b) Nucleon struction function in DIS. (c) Cross 
section for the nuclcon-nuclcon Drell-Yan process. 



2 Chiral-Odd Distributions hi^(x,Q 2 ) 

Inclusive hard processes cap be generally analyzed in the framework of the 
QCD factorization theorems. This theorem generalizes the idea of the Bjorken- 
Feynman's "parton model" and allows us to include QCD correction in a sys- 
tematic way. Here I restrict myself to the hard processes with the nucleon 
target, such as lepton- nucleon deep inelastic scattering (DIS, I + p — > I' + X), 
Drell-Yan (p + p' -> l + l~ + X), semi-inclusive DIS (I + p -> V + h + X). 
According to the above theorem, the cross section (or the nucleon structure 
function) for these processes can be factorized into a "soft part" and a "hard 
part" : The soft part represents the parton (quark or gluon) distribution in the 
nucleon and the hard part describes the short distance cross section between 
the parton and the external hard probe which is calculable within perturbation 
theory. For example, a nucleon structure function in DIS can be written as 
the imaginary part of the virtual photon-nucleon forward Compton scattering 
amplitude. (Fig. 1 (b)) According to the above theorem, in the Bjorken limit, 
i.e. Q 2 , v = P ■ q — ► oo with x — Q 2 /2v = finite, (Q 2 — —q 2 is the virtuality 
of the space- like photon, P is the nucleon's four momentum), the structure 
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function can be written as 

W(x,Q*) = J2 f ^-H^^,a s {^))^{y^ 2 ), (1) 
aJx y v m 

where $ represents a distribution of parton a in the nucleon and H a describes 
the short distance cross section of the parton a with the virtual photon, fj, 2 
is the factorization scale. In Fig. 1(b), <&" is identified by the dotted line. 
(Fig.l (a)). Similarly to DIS, the cross section for the nucleon- nucleon Drcll- 
Yan process can also be written in a factorized form at s = (Pa + Pb) 2 , Q 2 — > 
oo with a fixed Q 2 /s {Pa,b are the momenta of the two nucleons, Q is the 
momentum of the virtual photon): 

da ~£ f dy a fdy b H ab (—,—,Q 2 ;^,a s (fJ, 2 )) <J>"(y a , M 2 )<& 6 (y 6 , M 2 )(2) 

a b Jx a Jx b \ Va Vb M / 

where the two parton distributions, $ a and <E> b , for the beam and the target 
appear as was shown by dotted lines in Fig. 1(c). 

As is seen from Figs. 1(b), (c), the parton distribution can be regarded as a 
parton-nucleon forward scattering amplitude shown in Fig. 1 (a) which appear 
in several different hard processes. In particular, the quark distribution in the 
nucleon moving in the +63 direction can be written as the,light-cone Fourier 
transform of the quark correlation function in the nucleonfl 

r°° d?~ 

<t> a (x,v 2 ) = P + ^-e ixP - z (PS\r(0)rr(z)UPS), (3) 

where \PS) denotes the nucleon (mass M) state with momentum P^ and spin 
S^, and ip a is the quark field with flavor a. In (^), we have suppressed for 
simplicity the gauge link operator which ensures the gauge invariance and | M 
indicates the operator is renormalized at the scale /z 2 . A four vector a M is 
decomposed into two light-cone components a* = ^(a° ± a 3 ) and the trans- 
verse component a±. In (|h, z + = 0, z± = 0, and z 1 — 0. T generically 
represents 7-matrices, T — 7^, 7^75, cr^j,, 1. Q a (x : n 2 ) measures the distribu- 
tion of the parton a to carry the momentum k + = xP + in the nucleon, which 
is independent from particular hard proceeses. 

If one puts r = 7^,7^75, the chirality of $ and ip becomes the same, 
namely it defines the chiral-even distributions. Likewise, putting T = a^ v , 1 
defines the chiral-odd disributions. For the case of the deep-inelastic scattering 
(Fig. 1 (b)), the quark line emanating from the target nucleon comes back to 
the original nucleon after passing through the hard interactions. Since the 
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perturbative interaction in the standard model preserves the chirality except a 
tiny quark mass effect, the chirality of the two quark lines entering the nucleon 
in Fig. 1(b) is the same. Hence the DIS can probe only the chiral-even quark 
distributions. On the other hand, in the Drell-Yan process (Fig. 1 (c)), there 
is no correlation in chirality between two quark lines entering each nucleon. 
Therefore the Drell-Yan process probes both chiral even and odd distributions. 

The chiral-odd distributions h\(x, /i 2 ), h%(x, fi 2 ) in our interest are defined 
by putting T = <T^ij 5 in (|):t3 

/dz~ 
—e™ p -*(PS\r(0)a^ l5 r(z)UPS) 

= 2[h a 1 (x,^){S±^ Pv - S±„ Pli )/M 
+h a L (x, /j 2 )M(p fl n L , - p^n fl )(S ■ n) 
+h a 3 (x, fi 2 )M(S^n„ - S± u n ll )] (4) 

where we introduced two light-like vectors p, n (p 2 = n 2 = 0) by the relation 
P m = p M + Mi„M ; p _ n = l p - = n + = q if wc wr itc P+ = V,p= ^(1,0,0, 1), 

n = yf^(l; 0, 0, —1). V is a parameter which specifies the Lorentz frame of the 

system: V — > oo corresponds to the infinite momentum frame, and V — ► M/y/2 
the rest frame of the nucleon. is the transverse component of 5 M defined 
by S^ 1 = (S ■ n)p^ + (S ■ p)n^ + S^. One can show that $ a defined in (jjj) 
has a support — 1 < x < 1. If one replaces the quark field ip in (||) by its 
charge conjugation field C-ip T , it defines the anti-quark distribution $ a . In 
particular h® L 3 (x, /i 2 ) in (Q) are related to their anti-quark distribution by 

h t L ,3(- x ^ 2 ) = ~K,l,3( x ^ 2 )- 

$ a appears in a physical cross section in the form of the convolution with 

a short distance cross section in a parton level as is shown in ([!]) and (|^) . The 

cross section can be expanded in powers of "^t=j as 

M M 2 
a(Q 2 ) ~ ^(lnQ 2 ) + —B{\nQ 2 ) + -^C{\nQ 2 ) + • • • , (5) 

VQ 2 Q 

where each coefficient A, B, C receives logarithmic Q 2 -dependence due to 
the QCD radiative correction. In order to see how h\ l,3 can contribute in 
the expansion (|^), it is convenient to move into the infinite momentum frame 
(V ~ Q —* oo). In this limit the coefficient of h± t L,3 in (§|) behaves, respectively, 
as 0(Q), 0(1), 0(1/Q). Therefore if hi contributes to the A term in (0), 
can contribute at most to the B-term, and h,3 can contribute at most to the C- 
term. In general, when a distribution function contributes to hard processes at 
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spin 


average longitudinal transverse 


twist-2 
twist-3 


/i 9i hi 
e h L g T 



Table 1: Clasification of the quark distributions based on spin, twist and chirality. Under- 
lined distributions are chiral-odd. Others are chiral-even. 



most in the order of f } I , the distribution is called twist-r. Therefore 

h\, h,L, h$ in (0) is, respectively, twist-2, -3 and -4. 

Twist^2L distribution hi can be measured through the transversely polarized 
Drell-YanOi§Q0, semkinclusive deep inelastic-scatterings which detect pioro, 
polarized baryonsHErEl, correlated two pionsEx 

From the discussion above, one sees that it is generally difficult to isolate 
experimentally higher twist (r > 3) distributions in hard proceeses, since they 
are hidden by the leading twist-2 contribution (A term in (^J)). However, this is 
not the case for hL and gT- In particular spin asymmetries, they contribute to 
the B-term in±he absence of A-term: gx can be measured in the transversely 
polarized DISLL3, and hL appears in the longitudinal versus transverse spin 
asymmetry in the polarized nucleon-nucleon Drell-Yan processo. Therefore 
the Q 2 -evolution of <?t and Hl can be a new test of perturbativc QCD beyond 
the twist-2 level. 

Insertion of other 7-matrices in (||) defines other distributions, la Table 
1, we show the classification of the quark distributions up to twist-3 a There 
fi, gi,r, e is defined, respectively, by T — 7^,7^75,1 m.-(EJ). A similar clas- 
sification can also be extended to the gluon distributionstJT The distribution 
fi contributes to the spin averaged structure functions Fi^(x,Q 2 ) familiar in 
DIS. The helicity distribution gi contributes to the Gi(x,Q 2 ) structure func- 
tion measured in the longitudinally polarized DIS. By now there has been much 
accumulation of experimental data on fi and gi, and the data on ,01 triggered 
lots of theoretical discussion on the "origin of the nucleaa, spin" El. The first 
nonzero data on 92 (= .9t — gi) was also reported in Ref. tJ. 



3 Next-to-leading order (NLO) Q 2 -evolution of hi(x,Q 2 ) 

As we saw in the previous section, hi is the third and the final twist-2 quark 
distribution. It has a simple parton model interpretation as can be seen by the 
Fourier expansion of if) in (|j) . It measures the probability in the transversely 
polarized nucleon to find a quark polarized parallel to the nucleon spin minus 
the probability to find it oppositely polarized. Here the transverse polarization 
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' and has been frequently used for the 
The leading order (LO) Q 2 -evolution for hi has 



refers to the eigenstate of the transverse Pauli-Lubahski operator 75$* j_. If one 
replaces the transverse polarization by the longitudinal one, it becomes the 
helicity distribution g\. For nonrelativistic quarks, Mt j^) = gi(x,[i 2 ). A 
model calculation suggests, hi is the same order as gi.aEjO 

The Q 2 -evolution of hi is described by the usual DGLAP evolution equatioiO. 
Because of its chiral-odd nature it does not mix with gluon distributions. 
Therefore the Q 2 -dependence of hi is described by the same equation both 
for singlet and nonsinglet rlistrifyifrioris. For fi and gi, the NLO Q 2 evolu- 
tion was derived long time, 
analysis of experimentst 2 }! 

been known for some timeS. In the recent literature, the. pext-to- leading order 
(NLO) Q 2 -evolution has been completed by two papersE-a'Ej almost at the same 
time: Vogelsang^ presented the light-cone gauge calculation_£or thc-two-loop 
splitting function of hi in the formalism originally used for /p9. Wec3 carried 
out the Feynman gauge calculation of the two-loop anomalous dimension fol- 
lowing the method of Ref. Ila for fx. The results of these calculations in the 
MS scheme agreed completely. This result was subsequently confirmed by Ref. 
EO. In the following, I briefly discuss the characteristic feature of the NLO Q 2 
evolution of hi following Refs.E§E3. 

Analysis of (^) gives the connection between the n-th moment of hi and 
a tower of twist-2 operators: 



M n [hi(fi 2 )} 



dx x n hi(x, /J?) 



-1 
2M 



(PS ± \0»(fi 2 )S ± „\PS ± ) 



0» n = $a» a n a i l5 (in-D) n ^, 



(6) 



where S± stands for the transverse polarization and 0^(/i 2 ) indicates the op- 
erator O v n is renormalizaed at the scale /i 2 . The contraction with and Sj_ 
(recall Sj_ ■ n = 0, n 2 = 0) in (||) projects out the relevant twist-2 contribution 
from the composite operator. ("Twist" for local composite operators is defined 
as dimension minus spin.) By solving the renormalization group equation for 
0%S± V , one gets the NLO Q 2 dependence of M n [hi(fJ. 2 )] as 



M n [hi(Q 2 )} (a s (Q 2 ) 



M n [hi^ 2 )] W^ 2 ) 



a s (Q 2 ) - a s {n 2 ) Pi 




In \ 


^ fa 




" 2/3 ) 



where a s (Q 2 ) is the NLO QCD running coupling constant given by 



a s (Q 2 ) 

47T 



/3 ln(Q 2 /A 2 ) 



ftliiln(Q 2 /A 2 ) 
PoHQ 2 /A 2 ) 



(7) 



(8) 
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with the one-loop and two- loop coefficients of the j3- function fio = 11 — 2/3Nf 
and /3i = 102 — 38/3-ZV/ (Nf is the number of quark flavor) and the QCD scale 

parameter A. 7^ and jn are the one-loop and two-loop coefficients of the 
anomalous dimension j n for O^S^ defined as 



~ =^ 7 (o) + f^ N )%(D 
7,1 4tt 7 " + UJ ln 



(9) 



If one sets (3\ — > and 7A — > in (Q), the leading order (LO) Q 2 evolution 
is obtained, and jn are obtained, respectively, by calculating the one- 
loop and two-loop corrections to the two-point Green function which imbeds 
O v n S^v To obtain jn , calculation of 18 two-loop diagrams is required in the 
Feynman gauge. Since Ihe expression for jrP is quite complicated, we refer 
the readers to Refs. CJt3 for them. 



300 
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Figure 2: The IftiO anomalous dimension 7^' 1 ' in comparison with 7^ 9 ' 1 '. This figure is 
taken from Ref. c3. 



In order to get a rough idea about the NLO Q 2 dependence of hi, we 



plotted in Fig. 2 7 £ (1) (7^ for hi) 
nonsinglet fi and gi) for Nf — 3,5 



Ml) 

m comparison with 7„ 



(7^ for the 



One sees from Fig. 2 7^^ > jn 9 ^ 
especially at small n. This suggests that the NLO Q 2 evolution of hi is quite 
different from that of fi and gi in the small x region. The relation 7« > 
Jn 9 ^ is in parallel with and even more conspicuous than the LO anomalous 
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dimensions which read 

7* (0) = 2C F 

T M°) = 2C F 

To illustrate the generic feature of the Q 2 evolution, we have applied the ob- 
tained Q 2 evolution to a reference distribution for g\ and h\. As a reference 
distribution, we take GRSV g\ distributionEa and assume h\{x, fi 2 ) — gi(x, fi 2 ) 
at a low energy input scale (fi 2 = 0.23 GcV 2 for LQ-apd /i 2 = 0.34 GeV 2 for 
NLO evolution) as is suggested by a nucleon modelQ'lLa. We then evolve them 
to Q 2 = 20 GeV 2 and see how much deviation is produced between them. The 
result is shown in Fig. 3. As is expected from the anomalous dimension, the 
drastic difference in the Q 2 evolution between hi and gi is observed in thfi-s,mall 
x region, and this tendency is more significant for the NLO evolutionn3-E3L3 
(Although g\ for u-quark mixes with the gluon distribution, the same tendency 
in the difference from h\ is observed for the nonsinglet distribution.) 



n+l 



1 - 



(n+l)(n + 2) 



(10) 




Figure 3: (a) The LO Q 2 evolution of h± (denoted by Su) and gi (denoted by An) for the 
M-qupf|k. (b) The NLO Q 2 evolution of hi and g± for the u-quark. This figure is taken from 
Ref . m 



As another example, we showed in Fig. 4 the Q 2 evolution of the tensor 
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charge (Fig. 4(a)) and the first moment of hi and the nonsinglet /i (171) (Fig. 
4(b)). Although the NLO effect is sizable for the tensor charge, it is small for 
the first moment. (For the lattice QCD calculation of the tensor charge, see 
Ref. §.) 

In Ref. Ell, the Regge asymptotics of hi was studied and the small- a; be- 
havior was predicted to be hi(x) ~ constant (x — > 0). On the other hand, the 
rightmost singularity of jn and jn are, respectively, located at n = —2 
and n = — 1 in the complex n plane. Therefore inclusion of the NLO effect 
in the DGLAP asymptotics gives consistent behavior at x — > as the Regge 
asymptotics. This is in contrast to the (nonsinglet) fi and gi distributions, 
whose LO and NLO DGLAP asymptotics are the same. 




Q 2 ( GeV 2 ) Q 2 ( GeV 2 ) 

Figure 4: (a) The LO and NLO Q 2 evolution of the tensor charge normalized at Q 2 = 1 
GeV 2 . (b) The LO and NLO Q 2 evolu-Mpn of the first moment of hi and the nonsinglet /1 
(or gi). This figure is taken from Ref. £j. 

One of the interesting applications of the obtained NLO Q 2 evolution of 
hi is the preservation of the Soffer's inequalitypl 2\hf(x, Q 2 )\ < fi{x 7 Q 2 ) + 
g1(x-Q 2 ). Although the validity of this inequality hinges on schemes beyond 
LOeS, the NLO Q 2 evolution maintains the inequality at Q 2 > Q\ if it is 
satisfied at some (low) scale Q\ in suitably defined factorization schemes such 
as MS and Drell-Yan factorization schemes .cJ'cil. 

As was discussed in Sec. 2, a physical cross section is a convolution of a 
parton distribution and a short distance cross section. (See ([!]) and (||)) For the 
double transverse spin asymmetry (Att) in the Drell-Yan process, the NLO 
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short distance cross section has been calculated in Ref. El in the MS scheme. 
The analysis on Att combined with the NLO tranversity distribution predicts 
modest but not negligible NLO effect l£| 

4 Q 2 -evolution of h^(x,Q 2 ) and its N c — > oo limit 

In general, higher twist (r > 3) distributions represent quark-gluon correlation 
in the nucleon. Using the QCD equation of motion, one obtains from (0) the 
following relatioi 



h L (x,u 2 ) = 2x I ^hi(y,ti 2 )+h L (x,^), (ID 

Jx V 

~ iP+ r°° Ar~ r 1 r u 

h L (x,u 2 ) = — / ^ e - 2fap ' 2 / udu / tdi 

x (P5|| \^(uz)i l5 a tta gG l/ a (tz)z^z u ilj(-uz)\PS ll }, (12) 

where z 2 = 0, z + = and S\\ stands for the longitudinal polarization for 

the nucleon (5 P = Sff = — ^-n M ). This equation means that hi, consists 

of the twist-2 contribution and which represents quark-gluon correlation 
in the nucleon. We call the latter contribution "purely twist-3" contribution. 
(Expansion of ( |l2| ) produces twist-3 local operators, See ( p^[ ) below.) Equation 
( |Tl| ) reminds us of the Wandzura-Wilczek relationEEl for gT- 

pl dy_ 

V 



g T {x,(i 2 )= I —gi(y,n 2 )+g T (x,n 2 ). (13) 

Jx V 



For e and gT, one can write down relations similar to ( |l2|) . 

The Q 2 -evolution of the first and second terms in (|1 1|) is described sepa- 



rately. The evolution of Hl is quite complicated. A detailed analysis of (12) 
leads to the following relation for the n-th moment of hj&: 

[(n+D/2] 

M n [h L {n 2 )] = i-rXi OM^" 1 ^ ^^' (U) 

fc=2 ^ ' 

Rnk = \ $a Xa n x i l5 (in ■ D) k - 2 igG ua n» {in ■ D) n ~ k ^ - (k -f n - k + 2)] . 

(15) 

We note that the number of independent operators {R n k} (k = 2, • • • , [(n + 
l)/2]) increases with n. In the (Revolution, the mixing among {R n k\ oc- 
curs and the renormalization is described by the anomalous dimension matrix 
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[ln{g)\ki for {R n k}- If we put the LO anomalous dimension matrix for {R n k} 
as {"{n(g)]ki = {ct s /2Tr)[X n ]ki corresponding to (|^), the solution to the renor- 
malization group equation for {R n k\ takes the following matrix form: 



(PS|||# nfe (Q 2 )|PS|| 



[(n+l)/2] 

£ 

1=2 



A7 



(PSulRmi^lPSu 



(16) 



where L = ° 3 ^ 2 j . A„ for /i£ was derived in Ref. Ej. The Q 2 -evolution for 
TjT and e is also described .by matrix equation^imilar to. (^), and the solution 
was obtained in Refs.EH for gx and in Refo for e.CJ As is clear from (14) 
and ( |l6| ) M n [hL{Q 2 )] and A^ n [/ii(/Lt 2 )] are not connected by a simple equa- 
tion as in the case for the twist-2 distribution (see (@))li4 Although © gives 
complete prediction for the Q 2 evolution, it is generally difficult to distinguish 
contribution from many operators in the analysis of experiments. 
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Figure 5: (Right) Corrmlfite spectrum of the eigenvalues of the anomalous dimension matrix 
for h,L obtained in Ref. E£|. The symbol o denotes the one-loop anomalous dimension for hi. 
The solid line is the anomalous dimension (Ell) at large n. (Left) Spectrum of the eigenvalues 
of the anomalous dimension matrix for Hl at large Nc] The solid line denotes the analytic 
solution given in (|18[). This figure is taken from Ref. E3. 

In order to get some feeling on the Q 2 -evolution of hh : we plotted the 
eigenvalues of X n in Fig. 5 (right). For comparison, we also showed in the 
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same figure the LO anomalous dimension y„ /2 for hi- (Note the differene in 
convention between (f7|) and ([l6|).) As is clear from this figure, the Q 2 evolution 
of Hl is much faster than that of hi (See discussion below.) 

In the recent liter atureoEZlo, it has been shown that at large iV c (the 
number of colors), a great simplification occurs in the Q 2 - ev °hjtion of the 
twist-3 distributions. Recall X n in ( |l6| ) is a function of two Casimir operators 

C G = N c and C F = If one takes N c -> oo, i.e. C F -> N c /2, @ and 

(O) is reduced to 



.M„[M<2 2 )] = rS'OojWjX,.^ 2 )], (17) 




S-w.lEi-i + ^J. (is) 

This evolution equation is just like those for the twist- 2 distributions (see (0)). 
In Fig. 5 (left), we showed the distribution of the eigenvalues of X n obtained 
numerically at N c — > oo. The solid line is the analytic solution in (|l§|), which 
shows (|l8|) corresponds to the lowest eigenvalues at iV c — > oo. Since ( |l7| ) was 
obtained by a mere replacement — * N c /2 in (|l6|), the correction to the 
result is of 0(1/N 2 ) ~ 10 % level, which gives enough accuaracy for practical 
applications. 

This large- N c simplification is a consequence of the fact that the coefficients 
of R n k in (fTi) constitutes the left eigenvector of X n corresponding to the 
eigenvalue 7„ in this limit: 

[(«+i)/2] f , N / o; \ 

which implies that all the right eigenvectors of X n except the one corresponding 
to 7^ are orthogonal to the vector consisting of (l- This leads to ©. 

This large- ^^simplification of the Q 2 eiaplution was proved for the nons- 
inglet <?t in Ref.c3 and for hj, and e in RefO. The corresponding anomalous 
dimensions for gx and e are, respectively, 




(20) 



12 



Corresponding to three twist-3 distributions in table 1, there are three 
independent twist-3 fragmentation functionsxl (Their number is doubled-lo 6 if 
one includes final state interactions. See RefE) It has been shown in Refo that 
at large N c the Q 2 evolution of all these nonsinglet fragmentation functions 
is also described by a simple evolution equation similar to (|l7|). Therefore 
the simplification of the twist-3 evolution equation is universal to all twist-3 
nonsinglet distribution and fragmentation functions. 

To illustrate the actual Q 2 evolution of we have applied (|l7|) to the bag 
model calculation of /il. lZI (Fig. 6) Fig. 6(a) shows the bag calculation of /i Jj. 
At the bag scale, purely twist-3 contribution is comparable to the twist-2 
contribution. After the Q 2 evolution to Q 2 = 10 GeV 2 , Hl is dominated by the 
twist-2 contribution. This is a consequence of the large anomalous dimension 
( |l8| ) compared with the LO anomalous dimension of hiJj> in the right figure 
of Fig. 5). A similar calculation was done for gx in Refo. 



O 




1 .0 




Figure 6: (a) Bag model prediction for h^. The dashed line represents the twist-2 contribu- 
tion to hij. (b) Bag model prediction for evolved to Q?~r= 10 GeV 2 assuming the bag 
scale is fi 2 = 0.25 GeV 2 . These figures are taken from Ref. c3 

Another simplification of the twist-3 evolution occurs at n — > oo El In 
this limit, all the twist-3 distributions obey a simple DGLAP equation ( p~T[ ) 
with a common anomalous dimension which is slightly shifted from ( O ) and 



7„ = 4C F 




(21) 
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This evolution equation satisfies the complete evolution equation to the 0(ln(n) /n) 
accuracyEJ. In the right figure of Fig. 5, ( pl[ ) is shown by the solid line. One 
sees that it is close to the lowest eigenvalues except for small n. Combined 
with this n — > oo result, the large-iV c evolution equation in ( |l7| ) with (|l^) and 
( pp| ) for each distribution is valid to 0((l/N 2 )\n(n)/n) accuracy. 



5 Summary 

In this talk, I discussed the recent progress in perturbative QCD, especially the 
Q 2 evolution of the chiral-odd spin-dependent parton distributions h\{x,Q 2 ) 
and Hl{x, Q 2 ). 

The NLO Q 2 evolution for the transversity distribution h\(x, Q 2 ) was com- 
pleted in the MS scheme. This means the Q 2 evolution of all the twist-2 dis- 
tributions has been understood in the NLO level. The resulting Q 2 evolution 
of hi(x, Q 2 ) turned out to cause quite different behavior from the helicity dis- 
tribution 171(2;, Q 2 ) in the small x region if hi is assumed to be equal to g\ at 
a low energy scale. 

The LO Q 2 evolution for the twist-3 distribution hi,{x, Q 2 ) (and e(x, Q 2 )) 
was completed. Although their Q 2 evolution is quite complicated due to the 
mixing among increasing number of quark-gluon-quark operators, it obeys a 
simple DGLAP equation similar to the twist-2 distibution in the N c — > oo 
limit, as was the case for the Q 2 evolution of the nonsinglet gr(x, Q 2 ) distri- 
bution. The same simplification at N c — > oo was also proved for the twist-3 
fragmentation functions. Therefore this large- N c simplification was proved to 
be universal for the twist-3 distribution and fragmentation functions. 

By these studies, we now have at hand the necessary tools given by pertur- 
bative QCD for the analysis of the whole parton distributions. From the point 
of view of the "spin distributions" in the nucleon, one is more interested in 
the nonperturbative x dependence of those parton distributions. Armed with 
the development in the perturbative Q 2 dependence and the nonperturbative 
QCD techniques, we will be forced to challenge this issue when the new spin 
colliders start producing data on the new spin distributions. 
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